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FOREWORD 


The  two  computer  programs  described  in  this  report  were 
developed  by  Athena  Engineering  Company  under  Contract  DAAH01-76-C-0293 
with  the  U.S.  Army  Missile  Command,  Redstone  Arsenal,  Al^ama.  The 
development  was  sponsored  and  technically  monitored  by  the  Ground 
Equipment  and  Materials  (GEM)  Directorate,  US  Army  Missile  Research, 
Development,  and  Engineering  Laboratory.  Dr.  Bobby  Mullinix,  GEM 
Directorate,  served  as  the  Contracting  Officer's  Technical  Representative 
and  technically  monitored  all  work  on  the  project.  The  two  computer 
codes,  IMPLATE  and  IMPSHELL,  are  operational  on  the  MICOM  CDC  6600  com- 
puter system  and  on  a UNIVAC  1108  system  used  by  Athena  Engineering 
Company.  The  codes  are  written  in  Fortran  and  should  operate  satisfac- 
torily on  all  equivalent  computer  systems.  Copies  of  the  computer  code, 
either  on  cards  or  tape,  may  be  obtained  upon  request  from  the  contract 
sponsor. 
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1.  INTRODUCTION 


During  a nuclear  engagement  outside  the  atmosphere  a severe  threat 
to  the  structural  integrity  of  missile  structures  is  posed  by  the  x-ray 
component  of  the  radiant  energy  emission  from  a nuclear  weapon  detonation. 
When  the  levels  of  absorbed  energy  are  sufficiently  high,  sublimation  of 
the  outer  structural  surface  can  occur  with  a resulting  impulsive  load 
being  imparted  to  the  structure.  This  impulsive  load  can  produce  stress 
waves  in  the  structure  or  gross  structural  response  or  both.  Although 
the  influence  of  the  stress  waves  on  structural  integrity  has  been  studied 
widely  (1-3)  , the  influence  of  the  impulsive  loads  on  the  longer-time 
structural  response  has  received  less  attention  (4,5). 

The  objective  of  this  project  was  to  develop  an  ability  to  study  the 
structural  response  of  isotropic,  metallic,  flat  plate  and  cylindrical 
structures  subjected  to  impulsive  loads  produced  by  x-ray  induced  material 
sublimation.  It  is  assumed  in  this  study  that  the  material  remains  elastic 
and  that  the  x-rays  damage  the  structure  by  decreasing  the  structural 
thickness  where  sublimation  occurs.  Temperature  influence  on  material 
properties  is  neglected  for  simplicity. 

To  accomplish  this  objective,  a flat  plate  structural  response 
code,  QTRPLATE,  used  in  experimental  studies  at  the  U.S.  Army  Missile 
Command  (6)  was  modified  to  accommodate  x-ray  deposition  loads  and  a 
cylindrical  shell  code,  IMPSHELL,  developed  for  predicting  structural 
response  to  such  loads. 
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2.  IMPULSIVE  LOAD  CALCULATION 


The  radiant  energy  component  of  a nuclear  weapon  detonation  can 
sublimate  the  surface  layer  of  structures  exposed  to  the  energy  flux 
provided  the  fluence  is  sufficiently  high.  This  extremely  rapid  sublima- 
tion can  produce  impulsive  loads  on  the  structure  due  to  the  momentum 
transfer  resulting  from  the  rapid  movement  of  the  surface  particles. 

The  impulsive  loads  induce  structural  motion  and  associated  stresses, 
strains,  and  deformations. 


2. 1 General  Computational  Method 

The  Whitener  (11)  analytical  mode)  has  been  incorporated  into  botli  the 
flat  plate  and  circular  cylinder  structural  response  codes  described  in 
Sections  3.  and  4.,  respectively.  The  Whitener  model  requires  a knowledge 
of  the  energy  deposition  profile  through  the  thickness  of  the  exposed 
structure  along  with  appropriate  physical  properties.  The  model  can  be 
expressed  as 


where : 


Ax.  = 

3 


H . = 
3 

E = 
s 


n 


F = 
c 

(H-E 


* 2 

blowoff  impluse,  taps/cm 

3 

density  of  zone  j,  g/cm 

thickness  of  zone  j , cm 

total  energy  deposited  in  zone  j,  cal/g 

sublimation  energy  for  the  material,  cal/g 

number  of  zones 

conversion  factor,  4.186  X 10^  ergs/cal 
) . = excess  energy  in  zone  j 


it 


(2.1) 


A tap  is  1 bar  - microsecond. 


To  implement  this  model  the  distribution  of  the  deposited  energy  through 
the  thickness  of  the  structure  must  be  known. 
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2 . 2 KNISH  Computer  Code 

The  KNISH  photon  deposition  computer  program  (7)  used  to  compute  the  energy 
deposition  profiles  for  use  in  the  Whitener  model  has  evolved  from  a 
code  written  by  Dr.  John  Huntington  at  the  Moleculon  Company  for  the 
Defense  Atomic  Support  Agency  (DASA)  in  1968.  The  code  was  modified  by 
the  Physics  International  Company  in  1972.  The  code  treats  Compton  single 
scattering  and  photoelectric  energy  deposition  in  an  axisymmetric  slab 
geometry. 

The  input-output  operations  of  the  KNISH  code  have  been  modified 
for  use  in  the  codes  described  here  to  make  data  input  more  concise  and 
to  permit  direct  use  of  the  output  in  the  Whitener  model  for  impulse 
calculations.  The  following  simplifications  have  been  made  in  data  input: 


1)  the  spectrum  energy  mesh  is  computed  internally  by  a 
prescribed  formula  which  takes  into  account  the 
characteristic  temperature  of  the  nuclear  weapon 
detonation , 

2)  only  spectra  formed  by  superimposing  individual 
blackbody  spectra  are  permitted, 

3)  the  number  of  regions  in  the  structure  is  preset 
to  one , 

4)  the  number  of  dose  points  is  set  internally  in 
the  code  to  fifty. 


This  information  is  provided  to  relate  the  simplifications  made  in 
KNISH  input  procedures  to  the  original  code  described  in  Reference  7. 
Data  input  descriptions  for  the  flat  plate  code  IMPLATE  and  cylindrical 
shell  code  IMPSHELL  are  described  in  Sections  3.3  and  4.3,  respectively. 

The  accuracy  of  KNISH  was  investic-'‘-ed  by  comparing  the  accuracy 
of  energy  deposition  profiles  computed  witli  KNISH  with  standard  results. 
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Such  a comparison  for  copper  is  illustrated  in  Figure  2.1  for  various 
biackbody  temperatures.  As  can  be  seen  from  the  figure,  good  accuracy 
is  obtained  using  KNISH  for  all  biackbody  temperatures.  The  printed 
output  from  KNISH  for  this  case  is  contained  in  Table  2.1. 

2. 3 Flat  Plate  Model 

As  illustrated  in  Figure  2.2,  the  radiant  energy  striking  the 
flat  plate  is  assumed  to  occur  normal  to  the  plate  and  be  uniformly 
distributed.  For  purposes  of  impulse  calculation  the  plate  is  assumed 
divided  into  various  regions  or  zones  through  the  thickness.  These  zone 
descriptions  are  preset  in  KNISH  and  used  to  calculate  the  energy  deposi- 
tion profile.  From  a knowledge  of  this  profile  the  thickness  of  the 
sublimated  layer  is  determined  by  noting  where  in  the  structure  the 
deposited  energy  intensity  exceeds  the  sublimation  energy.  All  material 
absorbing  energy  at  levels  higher  than  the  sublimation  energy  will  be 
sublimated.  This  condition  is  illustrated  schematically  in  Figure  2.3. 

Once  the  deposition  profile  and  thickness  of  the  sublimated  layer 
is  determined,  the  Whitener  model,  Eq.  2.1,  can  be  used  to  compute  the 
total  impulse  in  taps  per  cm^.  The  impulse  is  assumed  imparted  to  the 
structure  by  a uniform  pressure  load  which  varies  linearly  with  time  as 
illustrated  in  Figure  2.4,  where  tg,  the  "shine  time",  is  input  as  data. 
The  maximum  pressure  P is  computed  by  use  of  the  equation 

P = 213/ts  (2.2) 

where  Ig  is  computed  using  Eq.  2.1. 

2. 4 Cylinder  Model 

Radiant  energy  is  assumed  to  impinge  on  the  cylindrical  shell 
model  as  indicated  in  Figure  2.5a.  Only  the  component  of  the  energy 
normal  to  the  surface  is  assumed  to  be  absorbed  in  the  cylindrical  struc- 
ture. This  causes  the  resulting  impulse  to  be  imparted  to  the  structure 
in  a nonuniform  manner  as  indicated  in  Figure  2.5b.  The  same  procedure 
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Figure  2.2. 


Schematic  Diagram  of  Impulsive 
Loading  of  a Flat  Plate 
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P(0,t) 

P(e,t) 


Radiant  Energy 


a)  Exposure  of  Shell  to  Radiant  Energy 


b)  Impulsive  Pressure  Distribution  on  Shell 
Figure  2.5.  Schematic  Loading  of  Cylindrical  Shell 
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used  in  the  flet  plate  model  to  calculate  the  impulse  is  used  for  the 
circular  cylindrical  model  except  the  nonuniform  impulse  distribution  is 
accounted  for  in  the  manner  indicated. 
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3.  FLAT  PLATE  RESPONSE  CODE 


3.1  Description 


The  computer  code  QTRPLATE,  developed  previous  to  this  work  (6), 
is  designed  to  predict  the  dynamic  response  of  flat  plates  to  air  blast 
loads  while,  simultaneously,  prestressed  by  inplane  thermal  and  mechanical 
loads.  The  code  implements  a finite  difference  solution  to  the  governing 
equations  of  motion  which  are  based  upon  small  deformation  theory  and 
linear  elastic  material  behavior.  A mode  superposition  and  shock  spectrum 
approach  is  used  to  compute  peak  values  for  the  stresses,  strains,  and 
displacements  in  the  plate.  Predicted  response  has  been  shown  to  agree 
closely  with  experimental  results  (8) . 

3. 2 Modifications 

To  convert  the  QTRPLATE  computer  code  to  accommodate  impulse  loads 
due  to  blowoff,  the  energy  deposition  code  KNISH  and  an  impulse  compu- 
tational routine  IMPULS  have  been  combined  with  QTRPLATE.  The  combination 
has  been  accomplished  by  creating  a separate  master  controlling  program, 
referred  to  as  IMPLATE,  which  simply  calls  the  subprograms  KNISH,  IMPULS, 
and  QTRPLT  to  perform  the  various  calculations  in  a sequential  manner. 

The  main  program  of  QTRPLATE  was  redesignated  as  QTRPLT.  The  integrity 
of  the  programs  KNISH  and  QTRPLATE  has  been  left  largely  undisturbed 
to  facilitate  ease  of  understanding.  Information  is  passed  between  the 
KNISH,  IMPULS  and  QTRPLT  subprograms  by  means  of  labeled  common  state- 
ments. An  organizational  chart  of  the  new  code  for  flat  plate  response 
IMPLATE  is  contained  in  Table  3.1. 

The  most  significant  modification  made  to  the  code  QTRPLATE  involved 
modification  of  the  subroutine  SPCTRM.  The  displacement  shock  spectrum, 
defined  for  each  gcjieralized  coordinate  of  the  flat  plate,  is  defined  as  (6) , 

..  . max  I 1 ^ 


i **^i  t>o  M . w.  / 
11' 

I o 


P(T)  sin  a)^(t-T)dT 


Table  3.1  ORGANIZATIONAL  CHART  FOR  IIQ’LATE 
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where  is  the  generalized  mass,  the  natural  frequency  associated  with 
the  i^.j^  degree  of  freedom,  and  P (t)  is  the  pressure  intensity.  Using 
the  expression  for  the  pressure  contained  in  Figure  2.4,  the  value  for 
(04)  can  be  expressed  as 


F.  (OJ.  ) 
1 1 


P (2+t  0).  ) 
o s 1 

t M.  0).  3 

s 1 1 


(2.4) 


The  value  of  FMAX  in  subroutine  SPCTRUM  has  been  modified  to  reflect  this 
change . 

In  addition  to  the  above  changes,  minor  modifications  have  been 
made  in  the  input/output  procedures.  These  changes  will  be  described  in  the 
following  section. 


3.3  Data  Input  Format 

Data  is  input  to  the  computer  code  IMPLATE  by  specifying  information 
on  punched  cards  as  follows: 

1.  Title  card 

2.  Fluence  card 

3.  Spectrum  definition  cards 

4.  Material  identification  card 

5.  Photoelectric  cross-section  cards 

6.  Plate  geometry  and  edge  restraint  card 

7.  Material  property  card 

8.  Mechanical  and  thermal  preload  card 

If  several  problems  are  to  be  run  sequentially,  caids  1 through  8 should 
be  prepared  for  each  problem  and  stacked  sequentially.  Normal  termination 
of  a run  stream  is  accomplished  by  placing  a card  with  the  letters  END 
in  the  first  three  columns  at  the  end  of  the  stacked  deck.  Details  of 
the  information  to  be  specified  in  cards  1 through  8 above  are  contained 
in  Table  3.2. 


13 


TABLE  3.2  DATA  CARD  DESCRIPTION  FOR  IMPLATE 


Card 

Columns 

Format 

Data  Item 

1 

1-80 

80A.1 

TITLE  - Any  alphanumeric  characters  chosen 

to  describe  problem 

2 

1-10 

ElO.O 

FLOONZ  - Fluence  (cal/cm^) 

11-12 

ElO.O 

SHINET  - "Shine  time"  (Ms) 

3 

1-10 

ElO.O 

TT  - Blackbody  temperature  (kev) 

11-12 

ElO.O 

FF  - Decimal  fraction  of  total 

fluence  at  temperature  TT 

Repeat  cards  of  type  3 for  as  many 
blackbody  spectra  needed  (5  maxium) . End 
cards  with  a blank  card. 

4 

1-10 

10A.1 

ZMAT  - Material  name 

11-20 

ElO.O 

DD  - Material  density  (g/cm^) 

21-30 

ElO.O 

ZZ  - Ratio  of  atomic  number  to 

atomic  weight 

31-40 

110 

NED  - No.  of  cards  containing 

photoelectric  cross-sections 

41-50 

ElO.O 

ES  - Sublimation  energy  (cal/gm) 

5 

1-10 

E12.0 

A1  - Coefficient  in  cross-section  equation 

11-20 

E12.0 

B1  - Coefficient  in  cross-section  equation 

21-30 

E12.0 

U - Upper  photon  energy  limit 

6 

1-10 

FIO.O 

A - Plate  width  (cm) 

11-20 

FIO.O 

H - Plate  thickness  (cm) 

21-26 

1A6 

BC  - Boundary  conditions  (CC  or  SS) 

7 

1-10 

FIO.O 

E - young's  modulus  (GN/m2) 

11-20 

FIO.O 

ANU  - Poisson's  ratio 

21-30 

FIO.O 

ALPHA  - Thermal  expansion  coefficient,  (l/^F) 

31-40 

FIO.O 

RHO  - Plate  density  (g/cm^) 

8 

1-10 

FIO.O 

TX  - X- tens ion  (N/cm) 

11-20 

FIO.O 

TY  - y-tension  (N/cm) 

21-30 

FIO.O 

TEMP  - temperature  (°F) 
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3.4  Sample  Problem 

To  illustrate  the  punched  card  format,  consider  the  problem  in 
which  a square,  aluminum  plate,  10  X 10  X 0.16  cm  is  subjected  to  an  inplane 
tension  preload  in  the  x-direction  of  magnitude  20  N/cm,  an  inplane 
tension  preload  in  the  y-direction  of  10  N/cm,  and  is  heated  to  a tempera- 
ture 20®F  above  the  reference  state.  Finally,  the  plate  is  subjected 

2 

to  a fluence  of  75  cal/cm  , from  a nuclear  weapon  with  a blackbody  spectrum 
and  a temperature  of  5 kev.  The  shine  time  of  the  nuclear  weapon  is  300. 
microseconds.  Young's  modulus  for  the  aluminum  is  68.95  GN/m2,  Poisson's 
ratio  is  0.3,  the  density  2.71  g/cm^,  and  the  thermal  expansion  coefficient 
is  4.0x10  ^/°F.  The  material  photoelectric  cross-section  and  atomic 
number  properties  are  contained  in  Table  3.3.  The  plate  is  assumed  to  have 
all  edges  clamped.  The  data  cards  to  run  this  problem  are  listed  in 
Table  3.3. 

The  computer  code  INFLATE  computes  the  natural  frequencies  and 
mode  shapes  for  the  prestressed  plate  and  uses  these  in  a spectral  analysis 
method  to  determine  a conservative  estimate  of  the  maximum  deflections, 
stresses,  and  strains  in  the  plate.  This  information  is  printed  out  along 
with  a description  of  the  input  data.  The  results  of  the  sample  problem 
described  above  are  presented  in  Table  3.4. 
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TABLE  3.3  DATA  CARDS  TO  EXECUTE  INFLATE  SAMPLE  PROBLEM 
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Table  3.4  (continued) 
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4.  CYLINDRICAL  SHELL  RESPONSE  CODE 


4.1  Description 

A computer  code,  designated  as  IMPSHELL,  has  been  developed  to 
predict  the  dynamic  response  of  a linear,  elastic  cylindrical  shell 
to  an  impulsive  load  produced  by  "blowoff".  The  cylindrical  shell  can 
be  subjected  to  an  internal  pressure.  A finite  difference  approximate 
solution  method  is  used  to  solve  the  governing  equations  of  motion. 

The  code  IMPSHELL  uses  a different  solution  technique  from  the  technique 
used  in  IMPLATE.  Whereas  the  mode  superposition  method  is  used  in  IMPLATE, 

a Newmark-Beta  (12)  time  intergration  procedure  is  used  in  IMPSHELL.  This 
permits  an  exact  determination  of  the  transient  response  of  the  shell  to 
be  computed-  for  comparison  with  experimental  data. 

4.2  Equations  of  Motion 

Applying  the  principle  of  virtual  displacements  to  the  shell 
structure  illustrated  in  Figure  4.1  permits  the  governing  equations 
of  motion  of  the  cylindrical  shell  in  the  v and  w directions  to  be 
expressed  as 


R(p-q)  (V-W,g)  + Ng  g ( R V,g  + W)  + 

+ Ng(v,gg  -f-  2w  g - V)  = phR^v  (4.1) 

and 

R(p-q)  (R-fv,g+w)  + NQ(w,gg-*!Sg-w-R) 

+ ^0'0  <'*''0"'')  = phR^w  (4.2) 
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q = external  pressure 


In  Eqs.  4.1  and  4.2,  p,  q,  v,  w,  R,  h,  and  0 are  defined  in  Figure  4.2. 

N is  the  force  per  unit  length  in  the  circumferential  direction,  and  p 
is  the  density.  Subscripts  preceded  by  a comma  denote  differentiation 
with  respect  to  a spatial  coordinate  and  dots  over  a variable  are  used  to 
designate  time  differentiation.  Eqs.  4.1  and  4.2  are  based  upon  the 
assumption  of  only  extensional  motion,  i.e.,  bending  effects  are  not 
included.  Such  effects  have  been  shown  by  Humphreys  and  Winter  (9) 
to  be  of  little  consequence  for  the  application  considered  here. 

A solution  to  Eqs.  4.1  and  4.2  has  been  obtained  by  linearizing 
the  equations  about  the  prestressed  state.  This  state  corresponds  to 
the  deformed  position  after  application  of  the  internal  pressure.  The 
resulting  linearized  equations  are  then  solved  using  standard  finite 
difference  techniques  for  both  the  space  and  time  variables.  A Nemark-Beta 
stepforward  solution  method  is  used  to  detcnnine  the  transient  response. 
The  organization  of  the  code  IMPSHELL  is  indicated  in  Table  4.1. 


4.3  Data  Input  Format 


Data  is  supplied  to  the  computer  code  IMPSHELL  by  specifying 
certain  information  on  the  following  punched  cards: 


1.  Title  card 

2.  Fluence  card 

3.  Spectriim  definition  cards 

4.  Material  identification  card 

5.  Photoelectric  cross-section  cards 

6.  Shell  geometry  and  internal  pressure  card 

7.  Structural  property  card 

8.  Time  and  print  sequence  card 


If  several  problems  are  to  be  run  sequentially,  cards  1 through  8 should 
be  prepared  for  each  problem  and  stacked  sequentially.  Normal  termination 
of  a run  stream  is  accomplished  by  placing  an  END  card  at  the  end  of  the 
stacked  data  deck.  Details  of  data  card  preparation  for  cards  1 through 
8 are  contained  in  Table  4.2. 
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Table  4.1  ORGANIZATIONAL  CHART  FOR  IMPSHELL 
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TABLE  4.2  DATA  CARD  DESCRIPTION  FOR  IMPSHELL 


Card 

Columns 

Format 

Data  I tom 

1 

1-80 

80A1 

TITLE  - Any  alphanumeric  characters  chosen 

to  describe  problem 

2 

1-10 

ElO.O 

FLOONZ  - Fluence  (cal/cm^) 

11-12 

ElO.O 

SHINET  - "Shine  time"  (ys) 

3 

1-10 

ElO.O 

TT  - Blackbody  temperature  (kev) 

11-12 

ElO.O 

FF  - Decimal  fraction  of  total 

fluence  at  temperature  TT 

Repeat  cards  of  type  3 for  as  many 

blackbody  spectra  needed  (5  maximum) . 

End  cards  3 with  a blank  card. 

4 

1-10 

lOAl 

ZMAT  - Material  name 

11-20 

ElO.O 

DD  - Material  density  (g/cm^) 

21-30 

ElO.O 

ZZ  - Ratio  of  Atomic  number  to 

atomic  weight 

31-40 

110 

NED  - No  of  cards  containing 

photoelectric  cross-sections 

41-50 

ElO.O 

ES  - Sublimation  energy  (cal/g) 

5 

1-10 

E12.0 

A1  - Coefficient  in  cross-section  equation 

11-20 

E12.0 

B1  - Coefficient  in  cross-section  equation 

21-30 

E12.0 

U - Upper  photon  energy  limit 

6 

1-5 

15 

NODES  - Number  of  nodes 

6-10 

15 

NOUT  - Increment  between  nodes 

selected  for  printing  results 

7 

1-10 

FIO.O 

DELT  - Time  increment  (ms) 

11-15 

15 

NINCS  - No.  of  time  steps 

16-20 

15 

DELP  - Time  steps  between  each  print  command 

8 

1-10 

FIO.O 

PINT  - Internal  pressure  (KN/m^) 
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Card 

Columns 

Format 

9 

1-15 

F15.0 

16-30 

F15.0 

31-45 

F15.0 

46-60 

F15.0 

61-75 

F15.0 

Data  Item 

R - Shell  radius  (cm) 

H - Shell  thic)cness  (cm) 

E - Young's  modulus  (GN/m^) 
ANU  - Poisson's  ratio 
RHO  - Density  (g/cm^) 
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4.4  Sample  Problem 


To  illustrate  the  data  preparation  steps,  consider  the  problem  of 

a cylindrical  aluminum  shell,  40  cm  in  radius,  0.6  cm  thick,  and  pressurized 

2 

internally  by  a pressure  of  6895  KN/m  subjected  to  the  energy  flux  from 

2 

a 5 kev  nuclear  weapon  detonation  such  that  the  fluence  is  75  cal/cm  . 

The  weapon  "shine  time"  is  300  microseconds.  Young's  modulus  for  the 

2 

aluminum  is  68.95  GN/m  , Poisson's  .ratio  is  0.3,  and  the  density 

3 

is  2.71  g/cm  . The  data  cards  required  to  execute  IMPSHELL  for  this 
problem  are  listed  in  Table  4.3. 

As  mentioned  previously,  IMPSHELL  computes  the  displacements, 
stresses,  and  strains  in  the  cylindrical  shell  as  a function  of  time. 

The  values  at  each  node  around  the  shell  are  printed  at  the  frequency 
specified  by  the  value  of  DELP.  The  results  of  the  sample  problem 
described  above  are  presented  in  Table  4.4. 
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Table  4.4  OUTPUT  FROM  IMPSHELL  SAMPLE  PROBLEM 


BEST  AVAIUBLLCQpy 


5.  CORRELATION  OF  EXPLODING  FOIL  IMPULSE  TESTING 
WITH  NUCLEAR  EFFECTS  LOADS 

The  impulsive  loads  produced  by  sublimation  effects  of  nuclear 
weapons  is  often  simulated  in  the  laboratory  by  the  use  of  "exploding 
foil"  techniques  (10) . In  such  test  programs,  thin  films  of  metallic 
foils  are  bonded  to  the  structural  surface  which  experiences  the 
impulsive  loads.  Large  doses  of  electrical  energy  are  discharged  into 
the  conducting  foils  from  capacitor  banks  resulting  in  sublimation  of 
the  foil  material.  This  sublimation,  in  turn,  loads  the  structure 
impulsively.  It  is  desirable  to  be  able  to  correlate  the  simulation 
loads  produced  by  the  exploding  foil  techniques  with  the  actual  loads 
produced  by  nuclear  weapon  effects. 

One  procedure  for  correlating  the  results  of  the  two  phencxnena  is 
based  on  the  equivalence  of  deposited  energy  and  total  impulse.  It  is 
assumed  that  the  source  of  energy  is  unimportant  and  only  the  amount  of 
energy  deposited  in  the  sublimated  layer  affects  the  impulse.  Action 
times  of  the  two  effects  are  assumed  to  be  approximately  the  same. 

The  procedure  for  calculating  the  impulsive  loads  due  to  nuclear 

weapon  effects  has  been  described  in  detail  in  Section  2.1.  This  procedure 

involves  first  determining  the  energy  deposition  profile  in  the  material 

due  to  exposure  to  radiant  energy  from  a nuclear  weapon.  This  profile 

is  calculated  using  an  energy  deposition  code  such  as  KNISH.  After  the 

profile  is  known,  the  total  impulse  is  calculated  using  the  Whitener 

model  described  in  Section  2.1.  A calculation  which  yields  the  total 

2 

impulse  per  unit  area  (taps/cm  ) produced  by  the  nuclear  weapon  effects. 

To  calculate  the  impulse  produced  by  an  exploding  foil,  use  is  made 
of  the  single-zone  approximate  method  for  impulse  calculation  defined  as 


where 
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2 

I = impulse  per  unit  area  of  foil  (taps/cm  ) 

® 3 

p = material  density  (g/cm  ) 

h = foil  thickness  (cm) 

E,  = energy  density  in  foil  from  capacitor  discharge  (cal/g) 
d 

E = sublimation  energy  of  material  (cal/g) 
s 


To  use  Eq.  5.1  to  calculate  the  impulse  produced  by  an  exploding  foil, 
assume  the  known  total  amount  of  energy  discharged  by  the  capacitor  bank, 
minus  any  losses,  is  deposited  uniformly  throughout  the  foil.  This 
allows  an  energy  density  E^  to  be  calculated  by  dividing  the  net  deposited 
energy  by  the  number  of  grains  of  foil  used.  Substitution  of  E^  into 
Eq.  5.1,  along  with  the  thickness  h and  the  material  properties  p and 
E^,  allows  the  impulse  per  unit  area  to  be  calculated. 

Conversely,  if  the  impulse  produced  by  a given  nuclear  effects 
loading  is  known,  the  energy  density  required  by  the  exploding  foil 
technique  to  produce  a similar  impulse  can  be  determined  by  solving 
Eq.  5. 1 for  E^,  i.e. , 


E^  = E + 
d s 


I 

B 

9150ph 


2 


(5.2) 


Eq.  5.2  has  been  evaluated  for  the  particular  case  of  aluminum 
(E^  = 3200  cal/gm,  p = 2.71  gm/cm^)  for  different  thicknesses  and  the 
results  presented  graphically  in  Figure  5.1.  This  illustration  permits 
a rapid  determination  to  be  made  of  the  energy  density  requirements  to 
simulate  the  results  for  a particular  nuclear  effects  produced  im- 
pulsive load. 

As  an  example  of  how  to  use  this  information,  assume  that  it  is 

desired  to  simulate  a nuclear  explosive  which  produced  an  impulsive  load 
2 

of  1000  taps/cm  . The  foil  material  to  be  sublimated  is  aluminum  with 
p = 2.71g/cm^  and  E^  =3200  cal/g  and  is  0.0508  mm  thick.  Evaluating 
Eq.  5.2  gives  a value  for  E^  of  3263  cal/g  required  to  sublimate 


Deposited  Energy,  log 


Impulse,  I (Kilo taps) 

Figure  5.1.  Correlation  of  Impulse  and  Deposited  Energy 


the  foil  at  the  desired  impulse  level.  If  the  structural  surface  area 

2 

to  be  loaded  impulsively  is  50  cm  , the  total  energy  requirement  is  ' 

= (3263)  (50)  (0.00508)  (2.71)  = 2,246  calories 

or,  = 9,401  joules.  Thus  9,401  joules  would  have  to  discharge  from 
the  capacitor  ban)c  into  the  foil  to  simulate  the  desired  result. 


6,  CLOSING  REMARKS 


The  computer  codes  described  in  this  report  were  designed  to 
meet  a particular  need  of  the  project  sponsor.  No  attempt  was  made  to 
create  general  purpose  computer  codes  capable  of  assessing  the  spectrum 
of  problems  associated  with  the  response  of  structures  to  electromagnetic 
radiation  from  a nuclear  weapon.  The  codes  serve  the  purpose  for  which 
they  were  intended:  correlating  specific  experimental  measurements  with 
specific  mathematical  models.  However,  due  to  limitations  placed  on  the 
models,  the  codes  are  not  presently  capable  of  assessing  the  vulnerability 
or  survivability  of  missle  structures.  Still  further  code  developments 
are  needed  to  address  these  problems.  The  codes  should  be  extended  to 
include  the  following  characteristics  to  make  them  more  accurate: 

1 . material  property  dependence  on  temperature 

At  the  high  temperatures  in  question  the  material  stiffness 
is  decreased  significantly  which,  in  turn,  significantly 
affects  the  structural  response. 

2.  large  deformation  effects 

The  loads  produced  by  nuclear  weapon  effects  .re  sufficiently 
severe  that  large  deformations  are  experienceu  before  the 
structure  fails. 

3.  failure  criteria 

Missile  structures  can  withstand  some  permanent  damage  before 
failure  occurs.  The  specific  details  of  structural  failure  is 
important. 

4 . elastic-plastic  material  behavior 

Material  yielding  and  plastic  deformations  occur  before  the 
structure  fails. 
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5. 


combined  blast,  impulsive,  and  operational  loads 
The  operational  loads  produced  by  acceleration  and  pressuriza- 
tion are  high  as  well  as  the  blast  and  impulsive  loads  produced 
by  nuclear  weapon  effects.  The  combined  effects  of  all  loads 
produces  failure. 

6.  arbitrary  geometrical  shapes 

More  complex  geometries  than  plates  and  cylinders  must  be 
evaluated. 

7.  composite  material  response 

The  response  of  layered,  reinforced,  composite  materials  to 
nuclear  effects  loads  is  important. 

Each  of  these  items  is  important  in  evaluating  missile  structural  integrity 
and  should  be  included  in  survivability/vulnerability  studies. 
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